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ABSTRACT 

Aims. The millisecond pulsar PSR J1740 — 5340 in the globular cluster NGC 6397 shows radio eclipses over ~ 40% of 
its binary orbit. A first Chandra observation revealed indications for the X-ray flux being orbit dependent as well. In 
this work, we analysed five datasets of archival Chandra data taken between 2000 and 2007 to investigate the emission 
across the pulsar's binary orbit. 

Methods. Utilizing archival Chandra observations of PSR J1740-5340, we performed a systematic timing and spectral 
analysis of this binary system. 

Results. Using a x 2 -test, the significance for intra-binary orbital modulation was found to be between 88.5% and 99.6%, 
depending on the number of phase bins used to construct the light curve. Applying the unbiased statistical Kolmogorov- 
Smirnov (KS) test did not indicate any significant intra-binary orbital modulation. However, comparing the counting 
rates observed at different epochs, a flux variability on times scales of days to years is indicated. The possible origin of 
the X-ray emission is discussed in a number of different scenarios. 
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1. Introduction 

The millisecond pulsar (MSP) PSR J1740-5340 (6397A), 
located in the core-collapsed globular cluster NGC 6397 at 
a distance of 2.5 kpc, was discovered during a systematic 
search for ms-pulsars in Galactic globular cluste rs using the 
Parkes Radio Telescope (jD'Amico et al.ll2001a|) . It is orbit- 
ing around a massive late type companion (> 0.14 M ; 
lOrosz fc van Kerkwiikll2003f) with an orbital period of 1.35 
day ( D'Amico et a l. 2001 q). The in clination of the binary 
syste m to the line of sight is i ~ 50° (|Orosz fc van Kerkwiikl 
|2003|) . The pulsar's spin period and period derivative are 
P = 3.65 ms (jD'Amico et al.l l2001al) and P = 4.0 x 
10" 20 s s _1 (jPossenti et al.ll2005f ). The contribution to P 
due to possible pulsar acceleration in the cluster's gravita- 



tional field is estimated to be smaller than 10 
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its large offset of 0.6 arcmin ( corresponding to eleve n core 
radii) from the cluster center (|D'Amico et al.ll2001bl ). The 
spin parameters imply a spin-down luminosity, a character- 
istic age and a dipole surface magnetic field of PSR J1740- 
5340 of E ~ 3.3 x 10 34 ergs s _1 , r c = P/2 P ~ 1.4 x 10 9 
years , and B± ~ 3.9 x 10 8 G, respectively (jPossenti et al.l 
120051) . 

The radio emission from PSR J1740-5340 is partially 
eclipsed in the orbital phase int erval 0.05 — 0.45 for ap- 
proximately 40% of each orbit (jD'Amico et alj l2001aHbh . 
Moreover, strong fluctuations of th e radio signals were ob - 
served at nearly all orbital phases (jD'Amico et al.|[2~001bl) . 
which led to the interpretation that the pulsar is orbit- 
ing within an extended envelope of matter released from 
the companion. The optical light curve of the companion 
showing tidal distortions provides strong evidence that PSR 



J 1740-5340 is orbiting a companion whose Roche lobe is 
nearly completely filled (jFerraro et al.|[200~lh . However, the 
system is in a radio-ejection phase in which accretion is 
inhibited by the radiation pressure exerted by the pulsar 
on the infalling matter. The strong interaction between the 
MSP flux and the plasma wind would explain the irreg- 
ularities seen in the radio signals from PSR J1740-5340 
(jFerraro et alj|200ll iBurderi et al.ll2002D . 

PSR J1740-5340 is in a wide orbit with a separation 
of ~ 6.5 Rq, so that the wind energy density impinging 
on the non-degenerate companion should be significantly 
less than that estimated for other much tighter eclipsing 
binary systems, such as the field ms-pulsar PSR B1957+20 
which has a separat i on of ~ 0.04 Rq from its companion 
(jArzoumanian et alj H99l . PSR J1740-5340 therefore is 
unlik ely to drive a wind of s ufficient density off its compan- 
ion (jD'Amico et al.ll2001bD . [O rosz fc van Kerkwiikl (2003) 
have studied the optical light curve of the companion and 
found no evidence of heating from the puls ar radiation, 
which supports the aforementioned inference. ISabbi et alj 
(2003) have investigated the chemical composition of the 
non-degenerate companion star and found a strong deple- 
tion of carbon. This suggests a scenario in which the com- 
panion is an evolved star that has lost most of its surface 
layers. In view o f the high rotation al velocity of the compan- 
ion (~ 50 km/s: ISabbi et al.ll2003l ). the stellar wind possibly 
can be strong enough to cause the mass-loss and result in 
an extended envelope of matter. 

X-ray emission from this binary syste m was detected 
with the Chandra X-ray Observatory (jGrindlav et al.l 
|2001| ) . Based on an observation in 2000 (cf. Table 1) which 
only covered ~ 40% of the binary orbit, iGrindlav et al.l 
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( 2002j) reported evidence that the count rate of the system 
appears to increase by a factor of ~ 2 at phase 0.4, just 
before the pulsar comes out of the radio eclipse. However, 
the short exposure time and the limited orbital phase cov- 
erage of this first observation (cf. top panel in Figure 1 and 
Table 1) did not support a detailed temporal and spectral 
analysis of the emission. 

In this publication we report on a search for X-ray or- 
bital modulation as well as on a spectral analysis of the PSR 
J 1740-5340 binary system making use of archival Chandra 
data that cover various orbital phases ranges. 

2. Observations and data analysis 

In total, five observations were targeted on the PSR J1740- 
5340 binary system using the Chandra Advanced CCD 
Imaging Spectrometer (ACIS). The first observation was 
performed on 2000 July 31 using the front-illuminated (FI) 
ACIS-I3 chip, while the other four observations taken on 
2002 May 13 and 15, 2007 June 26, and July 16 used the 
back illuminated (BI) chip ACIS-S3. We summarize the 
basic information of these observations in Table 1. Their 
binary orbit coverages are shown in Figure [TJ As can be 
seen from this figure, only Obs. IDs 7460 and 7461 cover a 
larger fraction of the binary orbit. In the search for intra- 
orbital flux modulation and spectral variation, we therefore 
considered only these two datasets which have higher pho- 
ton statistics and longer exposure times, while the other 
archival data were used to search for a possible flux vari- 
ability on time scales of days or years. Data analysis was 
restricted to the energy range 0.3 — 8.0 keV. Searching for 
X-ray pulses from PSR J1740-5340 was precluded by the 
inappropriate temporal resolution of the observing modes 
used. Standard processed level-2 data were used. Correction 
for aspect offset was applied before analysis. All data were 
processed with Chandra Interactive Analysis Observations 
(CIAO) version 3.4 software and CALDB version 3.4.2. 



2.1. Search for orbital modulation 

In order to search for a modulation of the X-ray flux as a 
function of orbital phase, we extracted the photons from 
a circle of 2 arcsecond radius centered on the radio pulsar 
position RA (J2000) = 17 ,l 40 m 44?59, Dec = -53 o 40'40'.'9 
(jD'Amico et al.ll2001bl ). In the 2007 observations, ~ 90% 
of all source counts are located within this cut radius. The 
photon arrival times were corrected to the solar system 
barycenter with the CIAO tool AXBARY (JPL DE200 
solar system ephemeris). W e used the pulsar ephemeris 
from iD'Amico et al.l (|2001bf ) and correct the photon ar- 
rival times for the orbital m otion of the binary system 
(|Blandford fc Teukolskvi [197(1 . Since PSR J1740-5340 is 
far out of the globular cluster core, the correction for grav- 
itational acceleration in the cluster potential was ignored. 

We tested the 2007 data for intra-orbital flux modu- 
lation by fitting a constant to the X-ray light curve. In 
order to have the signal-to-noise ratio for each phase bin 
higher than ~ 4, the numbers of phase bins were restricted 
to be within 10 to 15 bins per orbital period. In addition, 
we tested 20 and 40 bins per orbital period for compari- 
son. Using a % 2 -test, the significance for a flux modulation 
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Fig. 1. Orbital phase coverages of five Chandra observa- 
tions targeted on the binary system PSR J1740-5340. The 
Y-axis shows the X-ray source counts per bin per kilo- 
second in the energy range 0.3 — 8.0 keV. The X-axis is 
the bin ary orbital phase (from ephemeris of lD'Amico et all 
l2001bl ). </>=0.0 is at the ascending node of the pulsar orbit. 
Binning is so that the signal-to-noise ratio in each bin is 
> 3. Error bars represent the lcr uncertainties. All X-ray 
light curves are corrected for vignetting and the various 
instrumental sensitivities of the different observing modes. 
Labels are the Obs. IDs (cf. Table 1). 



over the observed orbit was found to be between 88.5% and 
99.6%, depending on the number of phase bins used to con- 
struct the light curve. Table [5] summarizes the results for 
various bin numbers. The X-ray light curve for 10 phase 
bins is shown in Figure [U 

The disadvantage of the \ 2 -test is its dependence 
on the number of phase bins. We therefore applied the 
Kolmogorov-Smirnov (KS) test to the unbinned light curve 
data in order to have a bin-independent statistical eval- 
uation of the X-ray emission variability (see Figure [3]). 
Calculating the corresponding KS probabilities between our 
data set and the cumulative distribution function generated 
by assuming a constant X-ray flux indicated no significant 
deviation between these two distributions. The significance 
for an intra-orbital flux modulation from this test is only at 
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Table 1. List of Chandra observations of PSR J1740-5340 in NGC 6397 



Obs. Date 


Obs. ID 


Instrument 


Data Mode 


Effective Exp. Time 


Orbital Phase 


Net count rate" 


Off-axis(') 


2000-07-31 


79 


ACIS-I 


Faint 


48.3 ks 


0.146 


-0.553 


1.89 ± 0.25 6 


0.8745 


2002-05-13 


2668 


ACIS-S 


Faint 


28.1 ks 


0.035 


- 0.267 


1.74 ±0.25 


0.9368 


2002-05-15 


2669 


ACIS-S 


Faint 


26.7 ks 


0.501 


-0.719 


2.89 ±0.33 


0.9214 


2007-06-22 


7461 


ACIS-S 


Very Faint 


90.0 ks 


0.189 


- 0.950 


2.26 ±0.16 


0.3732 


2007-07-16 


7460 


ACIS-S 


Very Faint 


149.6 ks 


0.437 


- 1.714 c 


2.92 ±0.14 


0.2129 



"ACIS-S net count rates in units of 10 3 counts/sec corrected to the on-axis. 

6 The ACIS-I on-axis net count rate is 1.31 ± 0.17. In order to derive a comparable ACIS-S net count rate from the observed 
ACIS-I rate we convert the ACIS-I rate by using the WebPIMMS tool (http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html). 
c The total exposure time covers more than one binary orbit. 



Table 2. Significance for intra-orbital flux modulation by 
fitting the light curves of different phase bins against a con- 
stant by using the x 2 test statistics. The best-fit mean levels 
are listed as well. 



Number of bins 


Mean level 


xi 


Signihcance(%) 


10 


37.1±2.7 


2.22 


98.2 


11 


33.4±2.4 


1.82 


94.8 


12 


29.6±2.1 


1.62 


91.2 


13 


28.6±2.1 


2.38 


99.6 


14 


26.4±2.0 


1.86 


97.0 


15 


24.7±1.9 


2.01 


98.7 


20 


18.4±1.1 


1.34 


89.0 


40 


9.2±0.5 


1.25 


88.5 
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Fig. 2. X-ray source counts (0.3-8.0 keV) vs. orbital bi- 
nary ph ase according to the ephemeris of iD'Amico et all 
(200iy). Two orbital cycles are shown for clarity. The back- 
ground noise level is found to be at ~ 0.35 counts/bin. 
4> = 0.0 corresponds to the ascending node of the pulsar 
orbit. Error bars indicate the la uncertainty. The shaded 
regions between the orbital phases 0.05-0.45 and 1.05-1.45 
mark the radio eclipse of the pulsar PSR J1740-5340. The 
dotted lines indicates the best-fit constant level and its la- 
uncertainty. 



the level of ^68%. A similar result is obtained if we restrict 
the analysis to the hard energy band above 2 keV. 



2.2. Spectral analysis 

To investigate whether the X-ray spectral behavior of PSR 
J 1740-5340 varies across the orbit, we analysed the X-ray 
spectrum within and outside the eclipsing region separately. 
The source and background spectra were extracted from the 
same region as mentioned in §2.1 and from a source-free re- 
gion close to the pulsar. Response files for the corresponding 
observations were created by using the tools MKAC1SRMF 
and MKARF of CI AO. The background-subtracted count 
rates from the PSR J1740 — 5340 binary system per ob- 
servation are given in Table [TJ Extracted spectra were dy- 
namically binned in accordance with the photon statistic 
in each dataset. Model spectra to Obs. IDs. 7460 and 7461 
datasets were fitted using XSPEC 11.3.2. There are 183±14 
and 457 ± 21 net counts in total for the spectral analysis 
inside and outside the eclipsing region, respectively. 

Assuming that the emission originates from the shock 
interaction of the pulsar wind with the wind of the com- 
panion star or from non-thermal emission processes in the 
pulsar magnetosphere, we expect the X-radiation to be syn- 
chrotron. To test this hypothesis, we fitted the spectrum 
with an absorbed power-law model (PL). For the PL fit, no 
significant variation of the spectral parameters is found be- 
tween the spectrum inside and outside the eclipsing region. 
A thermal bremsstrahlung model (TB) was also tested, 
which physically implies that the X-ray emission is either 
from a hot plasma presented in the binary system or from 
the free-free radiation of the companion's corona. We also 
tested whether a single blackbody (BB) can provide an ap- 
propriate modeling of the data. However, this BB model 
cannot provide any statistically acceptable description of 
the observed spectra (i.e. xt > 2) either inside the eclips- 
ing region or outside. 

The spectral parameters inferred from these fits are 
summarized in Table [3] The spectral parameters inferred 
from the PL /TB fits are found to be con sistent with those 
reported by iGrindlav et~aH (|200lL l2002f) . The best-fitted 
Nii for each spectrum varies marginally but is consistent 
with the values of ~ 1.03 and ~ 2.22 x 10 21 cm~ 2 inferred 
from the optical reddening of NGC 6397 and from the radio 
dispersion measure (DM) of PSR J1740-5340, respectively. 



3. Summary and discussion 

We have searched for the orbital modulation of the X- 
ray emission from PSR J1740-5340 associated with NGC 
6397. Only one of five archival datasets cover a full bi- 
nary orbit. Analysing this data with a Kolmogorov-Smirnov 
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Table 3. Parameters from spectral fits to the PSR J1740 — 5340 binary system. 



Orbital Phase 


Model a 


N H " 


T/kT(keV) 




F x c 




d 


Note 


0.05-0.45 


PL 
TB 


< 1.4 

< 0.8 


!- 8 -0.3 
7 n+10.4 
' - z -3.4 


1.39/20 
1.42/20 


1 Q+l-8 
± - a -0.7 
1 7+0.5 
1 - 1 -0.4 


1 4+1-3 
!-*-0.5 
1 q+0.4 
1 - J -0.3 


x 10 31 
x 10 31 


inside eclipsing region 


0.45-0.05 


PL 

TB 


1 S+ 1 - 

I O + 0.6 

L -°-0.9 


i 7 +o.i 


0.96/19 
1.00/19 


q q+1.4 
o -°-0.8 

2 9 +0 ' 4 


9 r+1.1 
z -°-0.6 
9 9+O.3 
ZiZ -0.4 


x 10 31 
x 10 31 


outside eclipsing region 


0.0-1.0 


PL 
TB 


6 +04 
< 0.7 


1 e+0.2 
l-°-0.1 

16.1± 23 6 2 


1.07/19 
1.11/19 


9 O+0.9 
z -°-0.8 

2 6 +0 ' 3 
z -°-0.4 


2 l +0 ' 7 
1 q+0.2 

l- y -0.3 


x 10 31 
x 10 31 


whole orbit 



Note: Col. 1 gives the orbital phase, Cols. 2-5 the spectral model, the hydrogen column density Nh, the best-fitted photon index 
or temperature, and the reduced x 2 together with degrees of freedom. Cols. 6-7 list the unabsorbed X-ray flux and luminosity 
in the 0.3 — 8.0 keV energy range. Quoted errors indicate the 68% confidence level for one parameter of interest. 

a. PL = powerlaw; TB = thermal bremsstrahlung 

b. In units of 10 21 cm -2 

c. Unabsorbed X-ray flux in units of 10 -14 ergs cm~ 2 s 

d. Unabsorbed X-ray luminosity in units of ergs s~ 
l|Grindlav et al.ll200ll i. 



and in the energy range of 0.3-8.0 keV. 
and in the energy range of 0.3-8.0 keV for a distance of 2.5 kpc 



Obs ID 7460 



0.1 
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Fig. 3. Kolmogorov-Smirnov test. The cumulative distri- 
bution function from the observed data is plotted in red 
and the blue dotted line shows the estimated cumulative 
distribution generated by assuming a constant X-ray flux. 

test revealed no significant intra-orbital flux modulation. 
However, correlating the ACIS-S vignetting-corrected net 
counting rates observed at various orbital phases in 2000, 
2002 and 2007 and assuming that the system in general 
shows no intra-orbital flux modulation reveals a ~ 3cr flux 
variability on time scales of days to years (cf. Col. 7 in Table 
1). Figure [4] depicts the ACIS-S net count rates observed in 
the various observations. 

In the following, we will discuss three scenarios which 
we believe could be the main source of the observed X-rays 
in the PSR J1740-5340 binary system. 



1, The Companion star 

The companion is observed to be a late-type star. It has an 
unusual position in the color-magnitude diagram (CMD) 



52500 53000 
MJD 



Fig. 4. ACIS-S3 count rates of the PSR J1740-5340 binary 
system vs. observation dates. 



as bein g as luminous as a main-sequence turn-off star bu t 
redder (jFerraro etail 120011 : lOrosz & van Kerkwiikl l2003l h 
This makes it quite similar to the "red straggler" ac- 
tive binaries which ar e identified to be X-ray sources in 
other globular clusters (lAlbrow et al.ll200H lEdmonds et al.l 
l2003aM iBassa et al.l [2008|). Detailed studies of X-ra; 



sour ces in various g lobul ar clusters by IBassa et al.l (|200 
and IVerbunt et all (|2008l ) have shown that the separatrix 
logL0.5-2.5kcV = 34.0 — 0.4My allows one to distinguish 
between the population of cataclysmic variables and that 
of magnetically active binaries. Here, My is the absolute 
V-band magnitude, which for the P SR J1740 — 5340 binary 
syste m was observed to be My ~ 5 (jOrosz fc van Kerkwiikl 
2003). Taking this together makes the properties observed 
from the companion star consistent with that from rapid 
rotating stars. The X-ray luminosities of late-type stars 
(i.e. F7 to M5) are found to be well correlated with their 
equatorial rotation al velocities Vmt , acco rding to Lx — 
10 27 u^ ot ergs s -1 (jPallavicini et all 11981 ). Ado pting the 
v rot ~ 50 km/s as reported bv lSabbi et al. (|2003l ) yields an 
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expected X-ray luminosity of ~ 2.5 x 10 30 ergs s _1 . Given 
that the Lx v s. v^t correlation has a scatter of one order 
of magnitude (jPallavicini et al.lll98ll ) it cannot be excluded 
that the X-ray emission observed from the PSR J1740-5340 
binary is coronal emission from the companion star only. 
No significant X-ray intra-orbital modulation is observed, 
which is unlike what we have seen in the optical bands. The 
optical emission is expected to come from the surface of the 
companion star while the X-ray emission comes from the 
corona. A s the companion star is inte rpreted to be tidal- 
distorted (|Orosz k, van Kerkwiikl [20031 ). the emitting area 
responsible for the optical intensity appears to be modu- 
lated. However, the X-ray emission is coronal and may be 
less or not affected by this distortion. 



Based on the preceding discussion, it is not possible 
to firmly distinguish the possible emission scenarios which 
may, alone or in combination, be responsible for the ob- 
served X-ray emission. To disentangle them would require 
detailed orbital-phase resolved spectroscopy along with a 
search for X-ray pulses from the millisecond pulsar. Further 
dedicated X-ray observations covering more binary orbits 
are certainly needed to further constrain the emission pro- 
cesses at work in this binary system. 

Acknowledgements . This work made use of the Chandra data archive. 
The first author thanks C.Y. Hui and M. Giidel for providing helpful 
suggestions and also acknowledges the receipt of funding provided by 
the Max-Planck Society in the frame work of the International Max- 
Planck Research School on Astrophysics at the University of Munich. 



2. The millisecond pulsar 

Considering a possible contribution from the pulsar mag- 
netosphere is justified by the pulsar's spin-down energy 
of E ~ 3.3 x 10 34 ergs s _1 . By assuming the best-fit X- 
ray conversion efficiency found for other field and millisec- 
ond pulsars, L x ~ 10~ 3 - 10~ 4 .E; (|Beckerll2009l) . the spin- 
down power of PSR J1740-5340 implies an X-ray luminos- 
ity which is of the order of ~ 10 30 — 10 31 ergs s _1 , in agree- 
ment with the observed luminosity (cf. Table 3) . Observing 
the pulsar in a mode suitable to resolve its X-ray pulses 
would be worthwhile and could constrain this interpreta- 
tion. 



3. Intrabinary shock emission 

iGrindlav et~aTl (|2002f ) suggested that the X -ray emission 
they detected from the PSR J1740-5340 binary system is 
due to the interaction between the pulsar and shocked gas 
lifted from the stellar companion. In the shock, the ex- 
tended matter in the orbit is compressed and gives rise 
to a power law distribution of synchrotron-emitting parti- 
cles, iV(7) cx 7 _p , by accelerating a fraction of these par- 
ticles to very high energies, where 7 is the Lorentz fac- 
tor of the emitting particles. This should lead to an en- 
hancement of X-ray emission from the system as the pul- 
sar enters the eclipsing region and a subseq uent decrease 
due to the increased absorption in this region. ICheng et al.l 
(2006D found that the X -ray luminosity of intrabinary 
shock emission in the 2-10 keV band can be described by 
L% ~ (0.3- 1.2) x 10 32 .D n (p - 1)/2 ^ +5)/4 ergs s"\ assum- 
ing a Lorentz factor for the synchrotron-emitting particles 
of 7 = (0.5 — 1.0) x 10 6 , the fractional energy density of 
electrons to be e e = 0.5 and the fractional energy density 
of the magnetic field to be es = 0.02 — 0.1. p and Dn de- 
note the index of the postshock electron energy distribution 
N("j) cx 7 _p , and the orbital separation in units of 10 11 cm. 
The X-ray photon index inferred from the spectrum is 
r ~ 1.8 which would suggest tha t the emission fro m the 
shock is in a slow cooling regime (ICheng et alj feooe). The 
photon index is related to p as T = (p+l)/2. The observed T 
implies that p is ~ 2.6. With £35 = 0.33 and D n = 4.5, Lf 
is found to be in the range of ~ (0.1 — 4.4) x 10 30 ergs s _1 , 
which is at the observed level. Intrabinary shock emission 
thus can explain the observed X-ray luminosity though the 
modulation depth and amplitude, which is also a function 
of the inclination angle, is unconstrained from this model. 
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